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A large scale extinction map of the Galactic Anticenter from 2MASS 
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ABSTRACT 

We present a 127° x 63° extinction map of the Anticenter of the Galaxy, based on 
(J — H) and {H — K) colour excess maps from 2MASS. This 8001 square degree map 
with a resolution of 4 arcminutes is provided as online material. The colour excess ratio 
{J — H) / {H — K) is used to determine the power law index of the reddening law (/3) for 
individual regions contained in the area (e.g. Orion, Perseus, Taurus, Auriga, Monoceros, 
Camelopardalis, Cassiopeia). On average we find a dominant value of 1.8 ±0.2 for the 
individual clouds, in agreement with the canonical value for the interstellar medium. We also 
show that there is an internal scatter of P values in these regions, and that in some areas more 
than one dominant (3 value is present. This indicates large scale variations in the dust prop- 
erties. The analysis of the Ay values within individual regions shows a change in the slope 
of the column density distribution with distance. This can either be attributed to a change in 
the governing physical processes in molecular clouds on spatial scales of about 1 pc or an Ay 
dilution with distance in our map. 

Key words: ISM: clouds - ISM: dust, extinction - Galaxy: structure- Infrared: ISM 



1 INTRODUCTION 

Understanding the formation of stars is inextricably linked to the 
formation, evolution and physical properties of molecular clouds. 
Dust is one of the best tracers of the distribution of material in giant 
molecular clouds. Its optical properties (reddening and dimming of 
star light) allow us to determine line of sight column densities of 
material. The knowledge of how clouds are structured (e.g. the dis- 
tribution of mass, size and column density of clumps within clouds) 
can be used to obtain constraints about how clouds fragment and 
will eventually lead to a better understan ding o f the initial mass 
function (e.g. Padoan, Nordlund, & Jones i \99f)). 

The mapping of the dust distribution in giant molecular clouds 
or entire cloud complexes becomes increasingly achievable due to 
the availability of all sky near infrared surveys such as 2MASS 
(Skrutskie et al. |200i)) as well as large computer clusters and Grid 
technology Basic techniques to map the colu mn de nsity of dust, 
i.e. extinction, a re: 1) star counting (e.g. Wolf (Il923h . Bok 0956), 
Froebrich et al. 2) colour excess (e.g. Lada et al. (11994 



Lomb ardi & Alves ( 200lh ): 3) or a combination thereof (e.g. Lom- 
bardi |200i); While most simple, the star counting techniques suf- 
fer from enhanced noise compared to the colour excess techniques, 
when applied to near infrared data. On the other hand, using colour 
excess methods to determine extinction requires the knowledge of 
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the reddening law, which might var y with position and/or column 
density (e.g. Froebrich & del Burgo (l2006l) ). In this pap er we a pply 
the most basic colour excess method from Lada et al. (Il994l) and 
determine the redd ening law in the near infrared following Froe- 
brich & del Burgo (l200d) . 



A large fraction of nearby giant molecular clouds is part of the 
Gould Belt. Towards the Anticenter of the Galaxy we find a variety 
of these clouds (e.g. Perseus, Taurus, Orion) in relative close prox- 
imity and with a relatively small amount of contamination from 
other clouds in the Galactic Plane. Additionally, a number of more 
distant cloud complexes (e.g. Cassiopeia, Camelopardalis, Auriga, 
Monoceros) are found in the same part of the sky We hence se- 
lected an area of 127° x 63° (in Galactic Coordinates) towards the 
Galactic Anticenter to be mapped in detail and to investigate the 
column density distribution in the major cloud complexes situated 
in this region. 



The paper is structured as follows. In Sect.|2]we describe our 
method to determine the extinction map. A description and discus- 
sion of the determination of the reddening law is given in Sect.[3] 
Finally we study the distribution of column densities in Sect.lH fol- 
lowed by the conclusions in Sect.|5] 
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Figure 1. Gray scale illustration of our 8001 square degrees extinction map of the Galactic Anticenter obtained from (J — H) and {H — K) colour excess 
determinations in 2MASS. Pixel values are scaled linearly from -0.33 mag to 10 mag optical extinction. The rectangles mark the regions selected for detailed 
analysis (numbers are identical to their ID in Table[T}, and magnified extinction maps of these areas can be found in Appendix lAll The full resolution FITS 
file of this map is available as online material. 



2 THE EXTINCTION MAP 

We selected all objects from the 2MASS point source catalogue 
with a Galactic Longitude of 116° < / < 243° and a Galactic 
Latitude of —32° < b < +31°. To ensure good photometric ac- 
curacy only objects with a quality flag better than 'E' are used. 
The area was then decomposed into 1° x 1° sized regions. Inde- 
pendent of the position we determined the median colour of stars 
(J—H and H—K) in 4' x4' sized boxes with an oversampling of 2. 
This leads to a pixel size of 2', 30x30 pixels per square degree, or 
3810x 1890 pixel for our 8001 square degree sized final map. The 
calculations required for such a project are embarrassingly parallel 
(i.e. the problem domain and the associated computational algo- 
rithms used to solve this problem can be optimally decomposed 
to use the available computer resources) and have been solved in 
a distributed fashion. We used the Grid-Ireland infrastructure (e.g. 
Coghlan et al. (l2005h ) to perform these calculations. 

The correct determination of extinction from colour excess 
methods requires that we remove the foreground stars to the cloud 
and that we reject young stellar objects. Both types of objects sys- 
tematically influence the determined extinction values. Foreground 
stars lead to on average bluer colours, hence the extinction is un- 
derestimated. Young stars are intrinsically red and hence lead to an 
overestimate of the extinction. Due to the large range of distances 
of the clouds in our map (ranging from 140 pc to at least 1 kpc) 
the removal of foreground stars is very difficult. In principle the 
number of foreground stars has to be evaluated for each individual 
cloud and removed statistically. The large number of clouds in the 
map with partly unknown distances, or overlapping clouds close to 
the Galactic Plane makes this an almost impossible task. However, 
if the stars in each box are not dominated by foreground objects and 
one uses the median colour of stars instead of the average, then the 
foreground stars are removed automatically. Froebrich & del Burgo 



(l2006h have shown (see their Fig. 9) that such an approach indeed 
reproduces the intrinsic extinction values of clouds extremely well, 
as long as the cloud is not too distant, or has a very high extinction. 
For very distant clouds or very high extinction regions, however, 
this approach detects zero extinction. Hence, locall y more sophis- 
ticated procedures such as descri bed in Lombardi (l2005h or sug- 
gested in Froebrich & del Burgo (l2006h have to be used to obtain 
correct extinction values for the highest extinction regions. Since 
only a very small portion of the map area is affected and we are 
investigating only the large scale distribution of clouds and their 
properties in this paper, the main conclusions of this paper are not 
significantly changed by this effect. However, when using the pre- 
sented extinction map one has to keep in mind that the extinction 
values in the very dense regions of clouds are underestimated and 
that there is a small Ay dilution for more distant clouds. Similarly 
to the removal of foreground stars, the use of the median colour also 
removes intrinsically red YSOs, as long as they do not dominate (as 
e.g. in the Orion Nebula cluster). 

The two median J — H and H — K colour maps of the entire 
field show small but significant large-scale gradients of the colour 
values of background (i.e. extinction free) regions. These gradients 
are caused by a change of the stellar population with distance from 
the Galactic Plane. In particular, the median colour of stars be- 
comes bluer at lower Galactic Latitudes and (less prominent) red- 
der towards Galactic Longitudes approaching the Anticenter of the 
Galaxy. To determine the (J — H) and {H — K) colour excess 
maps, we thus fitted a polynomial of 3rd order in / and b to the 
extinction free regions and subtracted it from the median colour 
maps. 

The resulting {J — H) and {H — K) colour excess maps are 
converted separately into H-band extinction maps using Equ. 8 and 
9 of Froebrich & del Burgo (l2006h . respectively. A power law index 
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Figure 2. Gray scale illustration the noise map. Pixel values show the one sigma noise level scaled linearly from 0.0 mag (white) to 0.6 mag (black) of optical 
extinction. The rectangles mark the regions selected for detailed analysis (numbers are identical to their ID in Table[T} 



of the reddening law of [3 = 1.73 is applied (see Sect.[3]for details). 
The resulting two H-band extinction maps are averaged and con- 
verted into our final optical extinctio n map by Av = 5.6S9 x Ah, 
using the reddening law from Mathis (ll99Qh . In Fig.[T]we show this 
map in gray-scale. Indicated in this figure are the regions of which 
we show magnified versions of the extinction map in Appendix lAll 
These are also the regions we analyse in detail in the forthcoming 
sections. The full resolution FITS version of our entire map will be 
provided as online material to this paper. Note that the construction 
of the final extinction maps as average of the Ah extinction maps 
from {J — H) and {H — K) colour excess maps leads to further 
dilution of the very high extinction cores. This will, however, effect 
the same areas that are dominated by foreground stars, as described 
above. Thus, the number of high column density regions in our 
analysis will be underestimated. 

The noise in the extinction map can in principle only be mea- 
sured in extinction free regions as the standard deviation of extinc- 
tion values from zero. To estimate the noise in all areas we have de- 
termined a noise map (see Fig.©. Since we used the median colour 
of stars in each pixel, we need to determine the standard error of 
the median for our noise map. For a large sample and normal dis- 
tribution this can in principle be done by cfmed = 1-25 • a/^/N, 
where a is the error of the colour of the individual stars and N the 
number of stars in each box. Since we use stars of different pho- 
tometric quality (see below) and we do not know if the sample of 
stars in each pixel follows a normal distribution, we cannot use this 
equation. However, the fraction of stars for any given photometric 
quality is very constant (see below). Thus, we can conclude that the 
noise in each pixel is proportional to 1/ VN for both, the (J — H) 
and {H — K) maps. We have thus scaled and co-added the 1 / ^/N 
maps of stars with (J and H) and (H and K) photometry, in order 
to obtain the final noise map. The scaling was done to match the 
measured noise in the Ay map in extinction free regions. As can 
be seen, the noise in the map varies significantly with Galactic Lati- 



tude, as well as with extinction. The 3 a noise ranges from 0.36 mag 
optical extinction, close to the Galactic Plane, to 1.0 mag Ay in the 
most northern and southern parts of the map as well as in high 
extinction regions. In Table[T]we list for all individual regions the 
average number of stars with detections in J and H, as well as the 
average noise in the area. There are pixels in the map (0.0066 %) 
where the noise is extremely large, i.e, there are no stars to deter- 
mine the colour. These regions are: 1) around very bright stars; 2) 
at |6| > 30°; 3) high extinction cores in Orion, Taurus, Perseus. 
For Ori A we determine an average noise of 0.24 mag, which can 
be compared to the 1 a noise level of about 0.2 mag Ay obtained 
by Lombardi & Alves (l200lh (see their Fig. 6) in this region using 
the optimised NICER method. 

There are further areas in the map where the extinction values 
are not very reliable. These are positions close to dense star clus- 
ters. Here, the high number of stars in the cluster dominates the 
extinction determination. Hence, depending on the colour of the 
cluster members we measure either a too low (mostly blue stars) or 
too high (mostly young, red stars) extinction. Since this is a com- 
mon effect, we over plotted circles at the positions of all known star 
clusters (found in SIMBAD) in the gray scale maps of our selected 
regions in Appendix lAll Additionally we over plot + signs at the 
positions of all possib le new star clusters with \b\ < 20° found by 
Froebrich et al. (l2007h . 

Does the use of 2MASS quality BCD stars instead of only 
quality A stars influence the noise in our map? Independent of the 
position about 40-45 % of all stars are of quality AAA. Given the 
larger errors in photometry for the quality BCD objects, one can 
estimate that the noise in the extinction map is about 10% (i.e. 
0.03 mag) higher, when including the lower quality objects . Ho\y - 
ever, as has been pointed out by Froebrich & del Burgo (l2006l) . 
systematic errors of extinction determinations based on colour ex- 
cess methods can be much higher than this, since the population of 
stars seen through a cloud can differ from the population of stars 
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in a control field (even if close-by). The more than twice as large 
number of stars per resolution element when including the lower 
quality objects will lead to a much better representation of the pop- 
ulation of stars and, as a result, to smaller systematic errors. This 
out weights in our opinion the slightly higher resulting statistical 
noise in our maps and justifies the use of the lower quality 2MASS 
sources. 

Note that the galaxies M31 (/=121.2, 6=-21.6) and M33 
(/=133.6, b=-313) show up as apparently high extinction regions. 
In M 3 1 the structure in our extinction map resembles quite closely 
the 170 /im emission map from Haas et al. (Il998h and in particular 
we are able to identify the dominant ring of material in the disc of 
the Andromeda Galaxy. 



3 THE REDDENING LAW 

The conversion of colour excess maps into extinction maps re- 
quires the knowledge of the extinction law. Our approach of using 
{J — H) and {H — K) colour excess maps enables us to deter- 
mine the power law index of the extinction law using the colour 
excess ratio R= {J — H)I{H — K) and 



R-l 



(1) 



from Froebrich & del Burgo (l2006h . We now can analyse the result- 
ing distribution of (3 values in our maps. For all the following anal- 
ysis of (3 we use only areas where both, the (J — H) and {H — K) 
colour excess values are 3 a above the noise. In Fig.[3]we show the 
histogram of the measured distribution of all [3 values (No{P)) as 
dotted line. It peaks at 1.73, in good agreenientwith the canonical 
value for the interstellar medium (Draine We hence use 

(3=1.1?) to convert our colour excess maps into extinction. The dis- 
tribution of (3 values shows a rather broad peak with a FWHM of 
1.6. This, seemingly, rather large value is in good agreement for 
the predicted scatter of about zb 40 % when using the colour ex- 
cess ra tio to determine [3 (see e.g. Fig. 1 in Froebrich & del Burgo 
(1200^) ). 

It is possible to investigate in more detail, if the measured dis- 
tribution of (3 values No{P) is in agreement with a single value of 
/3o=1.73 in the entire field or not. This can be done by convolving 
the observed distribution of signal to noise values, No{{J — H)), 
in our extinction map from Fig.[T]with the predictions for the width 
of the distribution of [3. We have, however, to allow for the possi- 
bility that there is an intrinsic scatter in the p values (c^^.^^) around 
(3o. W e hence have to adapt Equ. 1 1 from Froebrich & del Burgo 
(l2006h and obtain: 



2 



^(Ai-A2) 
(A1-A2)' 



1 + 

Q^co ^co V dR / 
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To determine the predicted distribution of (3 values we use: 



iVp(/3) = iV, 



{{J-H})-e ^{'fid-'-")) 



\{J- 



H) 
(3) 

were A^o is a normalisation factor and (3^ = 1.73 the peak of the 
measured (3 distribution. A Kolmogorov Smirnov test can then be 
performed with No(f3) and Np{f3) to determine the probability (P) 
that the two distributions are drawn from the same sample. 

In order to determine the predicted p distribution we require 



to estimate the parameters aco and 7co. The parameter aco rep- 
resents the square of the ratio of the noise in the two colour ex- 
cess maps and 7co represents the ratio of the square of the noise in 
one colour excess map and the covariance o f the tw o colour excess 
maps (see Equ.[2]and Froebrich & del Burgo (l2006h for details). For 
our dataset we measure values of: aco=3.6 and 7co=3.7. Note that 
the values for both parameters might slightly vary with position or 
hence from cloud to cloud. However, we have measured here the 
same values for aco and 7co as were obtain ed for 2MASS data in 
IC 1396 W by Froebrich & del Burgo (l2006l) . Due to this fact, and 
since the dependence of the predicted P distribution on aco and 7co 
is rather weak, we fix these values for all our determinations. 

This gives us the possibility to vary the intrinsic scatter cr^.^^ 
and the Po value until we find the highest probability P that the 
observed and predicted distributions of p values match. We list the 
obtained values of those parameters as well as the probabilities for 
the individual regions in Table[T] 

We over plot the resulting best predicted distribution for P val- 
ues of the entire field as solid line in Fig.[3] We find that with a 
98.9 % probability the two distributions are equal, when a value of 
po = 1.73 and an intrinsic scatter of cr|.^^ = 0.52 is assumed. 
The large intrinsic scatter shows, that there is a good chance that 
indeed variable p values are present in different parts of the field. 
Hence, repeating the same analysis for selected regions should lead 
to different P values in different clouds. If locally those clouds have 
more uniform dust properties, then a lower intrinsic scatter should 
be found. 

Thus, similar to the entire area, we investigate the peak value 
of the p distribution, the intrinsic scatter and how well the dis- 
tribution matches the assumption of a constant P value for our 
individual regions. In Table[T] we summarise the obtained results. 
We find an average value of p= 1.8 ±0.2. Almost in all investi- 
gated regions the P value is in agreement with the estimate for the 
entire field. Only in Auriga 2 (/3=2.05), Monoceros (P=l.52) and 
Camelopardalis 2 (P=\AO) significantly different values are found, 
hinting to different dust properties such as the grain size distribu- 
tion in these areas. Notably, regions with smaller p values tend to 
be closer to the Galactic Plane. However, it is not clear if this is 
caused by different dust properties or by the overlap of clouds in 
this area. Furthermore, small systematic off-sets in either (J — H) 
or {H — K) can cause changes in the determined p value. Thus, 
the procedure of fitting a polynomial to the extinction free regions 
in order to determine the colour excess maps might also contribute 
to these smaller p values. Since very close to the Galactic Plane 
there are no extinction free regions, the procedure might have in- 
troduced a small off-set in the colour excess maps and hence to the 
P values. Clearly, more detailed investigations are required to find 
the cause of the smaller P values very close to the Galactic plane. 

We also find that for most regions the intrinsic scatter of the P 
values is smaller or in the order of the value for the entire field. Only 
in Ori A is a significantly larger value found. The probabilities that 
measured and predicted values are drawn from the same sample are 
listed in Table[T] and the individual distributions are shown in Ap- 
pendix lA2l We find that in some regions the proposal of constant 
P values with the corresponding internal scatter is more probable 
than in others. Especially in the regions Monoceros and OriB low 
probabilities (P < 95 %) are found. We can conclude that in these 
regions we do not see a dominant fraction of the dust possessing 
constant optical properties. This might indicate large scale varia- 
tions of the optical dust properties within these clouds. 
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Table 1. Names and areas of our individual regions, average numbers of stars in the 4' x4' box, average noise in the field, peak in the distribution of f3, intrinsic 
scatter of the f3 distribution, probability that the measured distribution is in agreement with the predicted one for the assumption of a constant (3o value and 
intrinsic scatter in the whole region and slope of the column density distribution in the logN vs. plot. 



iU 


Name 


/ -range 


6-range 


=^ stars/box 




po 






slope 


1 


Auriga 1 


156-173 


-12 --03 


47 


0.15 


1.99 


0.45 


98.5 


-0.68 


2 


Auriga 2 


175-186 


-11 --02 


53 


0.14 


2.05 


0.39 


99.9 


-0.83 


3 


Cassiopeia 


116-127 


-04 -+10 


90 


0.11 


1.67 


0.43 


99.9 


-0.72 


4 


Camelopardalis 1 


127 - 148 


+03 -+15 


53 


0.15 


1.96 


0.49 


98.4 


-0.91 


5 


Camelopardalis 2 


137-152 


-05 - +04 


78 


0.12 


1.40 


0.54 


98.5 


-0.71 


6 


Camelopardalis 3 


152-161 


-04 - +08 


74 


0.12 


1.66 


0.49 


95.2 


-1.20 


7 


A-Ori 


188-201 


-18 --06 


36 


0.18 


1.87 


0.24 


96.3 


-0.81 


8 


Monoceros 


212-227 


-13 -+01 


67 


0.14 


1.52 


0.53 


82.7 


-0.68 


9 


Orion B 


201-211 


-17 --08 


34 


0.18 


1.77 


0.41 


91.2 


-0.38 


10 


Orion A 


208-219 


-21 --16 


20 


0.24 


1.66 


0.59 


99.8 


-0.28 


11 


Perseus 


156-163 


-25 --15 


18 


0.25 


1.97 


0.51 


99.9 


-0.30 


12 


Taurus 


164-178 


-19--10 


27 


0.20 


1.99 


0.47 


99.9 


-0.32 


13 


Taurusextended 


171-189 


-31 --17 


14 


0.29 


1.92 


0.43 


99.9 


-0.49 


14 


Entire Field 


116-243 


-32 -+31 






1.73 


0.52 


98.9 
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Figure 3. Histogram of the distribution of f3 values in the entire field (dot- 
ted line) measured for all pixels where the {J — H) and {H — K) colour 
excess values are at least 3 a above the noise. The solid line represents the 
predicted distribution for a value of [3=1.13, (7/3.^^^=0.52 and the mea- 
sured distribution of Ay values in the entire field (see text for details). A 
Kolmogorov Smirnov test shows that both distributions are identical with a 
probability of 98.9%. 

4 THE COLUMN DENSITY DISTRIBUTION 

In order to analyse and interpret the distribution of column density 
values in the entire map, as well as in the individual regions, we first 
have to recapitulate the selection effects induced by our method. As 
outlined in detail in Sect.|2l our method of determining the extinc- 
tion maps leads to three biases. (1) Very distant clouds, as well as 
very high extinction regions will not be detected by our method, 
since these areas are dominated by foreground stars. (2) There is 
a small Ay dilution with distance. (3) The averaging of extinction 
maps from {J — H) and {H — K) colour excess maps leads to di- 
lution of the high extinction cores. (4) Regions close to rich star 
clusters will be dominated by the colours of the cluster members. 
However, due to the use of the median stellar colours instead of 
the mean colours, all not too distant clouds without rich clusters 
and extremely dense cores, will have no systematically influenced 
extinction values. 

The general way of analysing the structure of molecular 
clouds is to determine the distribution of clumps within them. This 



includes the clump mass, size and column density distribution. We 
have hence used the clump find-ID algorithm (Williams et al. 
([1994) ) to extract all detectable clumps in our final extinction map 
and determined their size and mass distribution. Contrary to line or 
dust emission maps, extinction maps have real, very smooth large 
scale profiles which make it very difficult to extract structure infor- 
mation using an automated thresholding technique. This is reflected 
by the fact that the determined clump properties from our map de- 
pended, partly significantly, on the chosen levels and also on the 
pixel size of the map. 

Hence, instead of using the clump properties in the map, we 
will analyse the column density distribution in the individual re- 
gions. Since the distance of the clouds will influence the measured 
slope of the column density distribution (see below), we will not 
analyse the distribution of the entire map as a whole. We show in 
Appendix l A2 1 the Ay distributions for all individual regions (solid 
line). For each region we fit a linear polynomial to the log(N) vs. 
Av plots (dashed line) for extinction values larger than the 5 a 
noise (vertical line). The measured slopes are indicated in the in- 
dividual diagrams as well as in Table[T] 

The individual Ay distributions nicely show the above dis- 
cussed biases for high extinction values. Some show a clear 
drop-off for regions with very high column densities. Depend- 
ing on the distance of the cloud and its position, this drop- 
off occurs at different values. For OriA and OriB we see a 
clear drop at Ay = 10 mag, while for Perseus (Ay = 8 mag) and 
Camelopardalis 2 (Ay = 6 mag) lower values are found. The val- 
ues of the slopes range from about -0.3 to almost -1.0, with the 
exception of Camelopardalis 3 where -1.2 is found. 

The interpretation of these values is, however, rather difficult. 
Two effects could influence the measured slope: (1) the distance 
of the cloud; (2) the intrinsic distribution of material. The distance 
could influence the slope since in more distant clouds the extinc- 
tion per resolution element is averaged over a larger physical area 
in the cloud. Hence, for more distant clouds the column density 
distribution on larger physical scales is measured. One can artifi- 
cially move a cloud to larger distances by using a larger box- size 
for the extinction measurements. Indeed one finds that this results 
in a steeper column density distribution. This is in agreement with 
the values obtained from our data, where more distant clouds pos- 
sess on average steeper column density distributions. We basically 
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find that in all clouds where our resolution element is smaller than 
1 pc the slope of the distribution is -0.3. For all clouds were we are 
not able to resolve 1 pc (cloud distance larger than 850 pc) the slope 
in the column density distribution is about -0.75. This leads to the 
conclusion that there might be a change in the governing physics 
determining the column density distribution in molecular clouds on 
a scale of about 1 pc. Interestingly, this is about the Jeans length for 
a column density that is required for self- shielding and molecular 
hydrogen formation in molecular clouds (Hartmann et al. 42001) ). 
The already discussed small Ay dilution with distance could as 
well contribute to the change in the slope. Furthermore, there could 
also be intrinsic differences between clouds, as reflected by the dif- 
ferent slopes measured for e.g. Ori A and OriB, which are believed 
to be at the same distance. We note however, that if measured only 
for 6 < — 13° the slope for Ori B is -0.32, hence in very good agree- 
ment with Ori A and the other close-by clouds Taurus and Perseus. 
For the region with b > —13° in Ori B the slope is -0.75, indicating 
that these clouds are more distant. We further note that an overlap 
of different clouds, as may happen close to the Galactic plane, does 
not change the measured slope if the individual clouds are at about 
the same distance and have intrinsically the same slope. 



5 CONCLUSIONS 

Using Grid technology and 2MASS data we created the largest 
extinction map to date based on near infrared colour excess. Our 
127° X 63° sized map of the Galactic Anticenter contains the 
molecular cloud complexes of Orion, Taurus, Cepheus, Monoceros, 
Auriga, Camelopardalis and Cassiopeia. The three sigma noise 
level ranges from 0.36 mag optical extinction close to the Galactic 
Plane to 1.0 mag Ay at |6| ^ 30° and in high extinction regions. 
Our final optical extinction map has a resolution of 4' and a pixel 
size of 2'. 

The colour excess ratio (J — H)I{H — K) is used to deter- 
mine the power law index [3 of the wavelength dependence of the 
extinction in the near infrared. In the entire field the distribution of 
f3 peaks at a value of 1.73. The average value of the investigated 
cloud complexes is determined io [3 — 1.8 it 0.2. Both values are 
in agreement with the assumed values for the interstellar medium. 
We analyse in detail the distribution of [3 values in the entire field 
and individual regions. It is found that there is an internal scatter 
in the (3 values of about 0.4-0.5. A comparison of the measured 
and predicted (3 distribution by means of a Kolmogorov Smirnov 
tests shows that in most regions the proposal of a constant [3 with 
an internal scatter is valid. In some regions (e.g. Monoceros, OriB) 
low probabilities that observations and predictions match are found. 
Hence, there is no dominant (3 value in the cloud, a clear indication 
of large scale variations of the optical dust properties in these areas. 

We further analysed the column density distribution of the in- 
dividual regions. Due to our technique, the highest extinction cores 
{Av > 10 mag) cannot be identified reliably. For extinction values 
up to 10 mag Ay we find that the column density distributions show 
a linear behavior in the log(N) vs. Ay plot. On average a steeper 
column density distribution for more distant clouds is found, i.e. the 
Ay distribution is steeper when measured on larger spatial scales. 
This seems to indicate a change in the governing physics of the col- 
umn density at scales of about 1 pc or an Ay dilution with distance 
in our map. 
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APPENDIX A: INDIVIDUAL REGIONS 
Al Extinction and (3 maps 

In the following we show for all individual regions investigated in 
this paper larger versions of the extinction and 13 maps. Note that 
the full resolution FITS file of the extinction map is provided as 
online material to the paper. We overplot in the maps the position 
of all known star clusters in the field found in SIMBAD (circles) as 
well a s the n ew star cluster candidates (+ signs) found in Froebrich 
et al. (l2007h . in order to mark regions were potentially the cluster 
members might dominate the extinction determination (see Sect.[2l 
for details). 
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Figure Al. Left panel: Gray scale extinction of the region Auriga 1. Extinction values are scaled linearly from -0.33 to 7 mag of optical extinction. Right 
panel: Gray scale map of [3 values in the region Auriga 1 for areas where the < J — if > and < H — K > values are more than 3a above the noise (for the 
remaining areas a value of ^5=1. 73 is adopted). The [3 values are scaled linearly from zero (black) to three (white). In both panels we indicate the positions of 
known star clusters from the 2MASS database as circles and from the list from Froebrich et al. (2006) as crosses. 
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Figure A2. As Fig. lA2l but for the region Auriga 2. Extinction values are scaled linearly from -0.33 to 6 mag of optical extinction. 
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Figure A3. As Fig. lA2l but for the region Cassiopeia. Extinction values are scaled linearly from -0.33 to 7 mag of optical extinction. 
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Figure A4. As Fig. lA2l but for the region Camelopardalis 1. Extinction values are scaled linearly from -0.33 to 7 mag of optical extinction. 
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Figure A5. As Fig. lA2l but for the region Camelopardalis 2. Extinction values are scaled linearly from -0.33 to 8 mag of optical extinction. 
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Figure A8. As Fig.|A2lbut for the region Monoceros. Extinction values are scaled linearly from -0.33 to 7 mag of optical extinction. 
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Figure A9. As Fig. lA2l but for the region OriB. Extinction values are scaled linearly from -0.33 to 10 mag of optical extinction. 
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Figure AlO. As Fig.|A2]but for the region Ori A. Extinction values are scaled linearly from -0.33 to 12 mag of optical extinction. 
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Figure A13. As Fig. lA2l but for the region Taurusextended. Extinction values are scaled linearly from -0.33 to 5 mag of optical extinction. 



A large scale extinction map of the Galactic Anticenter from 
A2 f3 and column density distributions 

In the following we show for all individual regions investigated in 
this paper the distribution of p and column density values. 
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Figure A14. Left: Measured column density distribution in the region Auriga 1 (solid line). As dashed lined we overplot the linear fit to the area above the 
5 (7 noise level (vertical line). The slope of this fit is indicated in the graph. Right: Measured distribution of /3 values in the region Auriga 1 (dotted line) and 
overplotted model prediction for the assumption of a single ^ value with internal scatter in the entire region. The measured distribution includes only pixels 
where the < J — if > and < H — K > values are more than 3a above the noise. 
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Figure A15. As Fig. lA14l but for region Auriga 2. 





Figure A16. 



As Fig. lAlU but for region Cassiopeia. 
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Figure A19. As Fig. lA14l but for region Camelopardalis 3. 
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Figure A22. 



As Fig. lAlD but for region OriB. 
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Figure A23. As Fig. lAT4l but for region Ori A. 
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Figure A24. As Fig. lA14l but for region Perseus. 





Figure A25. As Fig. lA14l but for region Taurus. 
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Figure A26. As Fig. lA14l but for region Taurusextended. 



